Generalized waveguide approach to tight-binding wires: Understanding large vortex 

currents in quantum rings 
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We generalize the quantum waveguide approach to Hiickel or tight-binding models relevant to 
unsaturated tt molecular devices. A Landauer-like formula for the current density through internal 
bonds is also derived which allows for defining a local conductance. The approach is employed to 
study internal circular currents in two-terminal rings. We show how to predict the occurrence and 
the nature of large vortex currents in coincidence with vanishingly small currents in the leads. We 
also prove a remarkably simple formula for the onset of a vortex regime. 
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Tlie Quantum Waveguide Approach (WGA) intro- 
duced in Refs. [l|H has been extensively used to study 
multiterminal mesoscopic structures, like quantum rings 
or quantum wires, possibly in the presence of impuri- 
ties, magnetic fields, Rashba or Dresselhaus interactions. 
The basic idea of the WGA is to calculate the one- 
particle wavefunction ipn{x) 



ikx 



-ikx 



in the 

n-th wire by imposing the continuity at each vertex, i.e., 
</3i(0) = (^2(0) = . . . = ¥'A^(0) {M being the number of in- 
tersecting wires at the vertex), and the additional condi- 
tion Vn(0) = 0, which implies current conservation 
at the vertex (the derivative is taken along the incoming 
direction). It is worth noticing that current conservation 
is actually fulfilled by the weaker condition 



E^»(o) 



-<Pi(0), 



(1) 



where r is an arbitrary real number. To our knowledge 
such arbitrariness has never been discussed. 

In recent years it became possible to attach aromatic 
molecules or atomic chains to leads. These structures 
arc geometrically similar to their mesoscopic counter- 
parts and call for a Hiickel-like or tight-binding (TB) 
generalization of the WGA due to the inadequacy of 
the continuum free-particle description. In the discrete 
case the difficulty stems from the fact that we cannot 
impose a condition on the derivative of the wavefunc- 
tion. Approaches based on Green's functions^, iterative 
procedures^ or source-and-sink-potcntials^ have been 
proposed but none of them is directly related to the 
WGA. 

It is the purpose of this work to show how to general- 
ize the WGA to the TB case (TBWGA) and to use the 
method to predict the occurrence of large vortex currents 
in quantum rings observed for the first time in Ref. 3- 
The continuum case is recovered by a proper limit of the 
TB parameters and allows us to understand the physical 
meaning of the real constant r in Eq. ([1]). 




FIG. 1: Illustration of a vertex with A/" = 4. 



Tight binding WGA: We consider a generic system con- 
sisting of TB chains with at least one end-point in com- 
mon. In Fig. [T]wc illustrate a vertex of the system with 
TV = 4 intersecting chains. Let e„ be the onsite energy of 
the Ti-th chain, t„ the hopping parameter between near- 
est neighbor sites and ey the onsite energy of the vertex. 
We denote with ipenU) = flenC*'''"^ -I- 6e„e~*'^"-' the ampli- 
tude on the j-th site of the n-th chain of an eigenstate of 
energy e ~ en + 2t„ cos fc„. As in the continuum case the 
continuity of the wave function yields Af —1 independent 
equations 



V'el(O) = Ve2(0) 



(2) 



The additional condition, which plays the role of Eq. ([T]) 
in the WGA, is obtained by projecting on the vertex site 
the stationary Schrodinger equation: 



,(l) + £vVsl(0) =eVei(0). 



(3) 



Equations ([2|) and ([3|) provide TV independent equations 
for each vertex of the system. Thus, for a system having 
V vertices with TV, i = 1, . . . ,V, intersecting chains at the 
i-th vertex the above procedure yields '^^J^i equations. 
Letting V be the number of chains with both end-points 
belonging to the set of vertices and Q be the number 
of semiinfinite chains with one end-point connected to a 
vertex, i.e., the number of terminals, we have 2V + Q = 
'Ylii^i s-nd hence a degeneracy D < Q for each energy 



2 



level (the number of unknown constants {a^n}-, {ben} is 
27^ + Q + D where D is the number of terminals for which 
|£ — £„| < 2|t„|), as it should be. 

To recover the WGA we employ a three-point dis- 
cretization of the kinetic term (our argument does not 
rely on this specific way of discretizing). Then, i„ = 
— 1/(2A^) and e„ = 1/A^, where A is the spacing be- 
tween two points of the continuum wire (we use atomic 
units). The amplitudes i/'eii(l) correspond to the am- 
plitudes of the wavefunction at a distance A from the 
vertex. For clarity we then rename ipenU) with ip^nU'^) 
and rewrite Eq. ^ as 



Af r 



E 



V'en(A) -1^e«(0) 



A£>Vel(0) - -2eA7/.el(0), 



(4) 

where ey = 2£y — A/'/A^. Taking the continuum limit 
A — > we recover Eq. ([T]) provided that liniA^o A£y ~ 
r, which implies sy ^ r/A. Thus, the constant r is the 
amplitude of the J-like potential rS{x) at the vertex and 
is zero only for smooth potentials. 

Two terminal systems: We now focus on two termi- 
nal systems and obtain a Landauer-like formula for the 
current density through a generic bond. Let fl"(°) = 
-ff leads + Hdev + ^^tun bc the Hamiltonian of the system in 
equilibrium. The Hamiltonian of the left (L) and right 
(R) leads is 



Hi 



leads 



t 



3<0 



(cjcj^i +h.c.) +t V(4c,+i +h.c.), (5) 



j>0 



with nearest neighbor hopping t. The device is de- 
scribed by Hdev = Y.nm=i tnmdidm with real parame- 
ters tnm = tmn and TV the total number of sites. The 
device is connected to the left lead through site 1 and to 
the right lead through site M < N (see, e.g.. Fig. [5]). 
The tunneling Hamiltonian is 



tL{d\c^i + h.c.) + tnidljci + h.c). (6) 



We are interested in the long-time limit of the current 
density when an external bias Ua is imposed on lead a ~ 
L,R. 

At zero temperature the long-time limit of the density 
matrix Px,x', with x,x' site indices of either the leads or 
the device, is given by the sum of a steady-state contri- 
bution p^fl, and a dynamical contribution^ The steady- 
state contribution can be written in terms of left-going 
eigenstates IV'ei?) and right-going eigenstates l^/'ei) as 



(s) 

Px.x' 



E 



2|t| + (7c 



de 
2^ 



^e^{x')rec.{x). (7) 



with £f the equilibrium Fermi energy. The states are 
normalized according to {ipea\i^e' p) = 2TrSaf3S{e — e'). 
Without loss of generality we choose Ul > Uji and split 
Eq. ([7]) into three terms containing the contribution of 
the left-going evanescent states with energy in the range 



{—2\t\ + Uii,—2\t\ + UL); the left- and right-going current- 
carrying states with energy in the range (— 2|t| + UL,eF + 
Ur), and the right-going current-carrying states with 
energy in the range (£f + Ur.Ey + Ul)- The evanes- 
cent states can be chosen real-valued since they are non- 
degenerate and the biased Hamiltonian is invariant under 

is) 

time-reversal. Thus, the imaginary part of ^, simplifies 
to 



x.x' 



^[i^.^ix'Wed^)] 

-Y^^\i,Ux'WM. (8) 
-2|tn-c/z. 27r ^ 



Let us first consider the case in which both x, x' are site 
indices of the same lead. The amplitude on the leads 
tpeaij) = {0\cj\ipea) of a normalized scattering state is 



ReLe~ 



TeLe"i' 



j <0 

j <0 
J >0 



, (9) 
,(10) 



with £ = 2tcos{q) + Ul = 2tcos{q) + Ur, and Vais) = 
l/-\/4t2 — (s ~ UaY ttic density of states in lead a. Ex- 
ploiting current conservation and the orthogonality con- 
dition between left- and right-going eigenstates, i.e., 
R*eL^eR/ ^hii^) + RsrT*!^/ vr{£) = 0, it is straightforward 
to show that for a; = j > and x' = j ^ 1 the second 
term in Eq. ([8]) vanishes while the first term reduces to 
the well known Landauer formula 



Me) 
vr{£) ' 



(11) 



Below we show that the possibility of expressing ^[p^'^^/] 
as an integral over the bias window is valid for all sites 
including those in the central device. Let us express i/'efl 
as a linear combination of ip^a and the time-reversal state 

ib'^ =th* 



i'eRix) 




^T^^ [^eL(^) - KlAUx)] . (12) 



Extracting the transmission and reflection coefficients 
T^R and Ri;r one can easily verify that ip^Rix) is or- 
thogonal to ^Psl{x). Inserting Eq. (fT2|) into Eq. 
([5]) and exploiting current conservation for the right- 
going scattering state, one realizes that the imaginary 
part of '^a'ipeaix')'4'*aix) is identically zero and hence 

only states in the bias window contribute to 3[p^'^2']- 
Thus, the long-time limit of the current density Jnm = 
2tnm'^[pn,m] through an internal bond of the device con- 
necting site n to site m can be expressed in a Landauer- 
like formula. In linear response Jnm = GnmiUh — Ur) 
and exploiting the above result the local conductance 
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G is given by 

Gnm = "^""l. = —Q[%peFL{mW^^L{n)]. (13) 

We next specialize the analysis to devices consisting of a 
TB ring and address the existence of vortex regimes. 




FIG. 2: Illustration of a ring device with M — 5 and N = 8. 
The last site of the lower arm is M' = N — M. 

Conductances in two terminal rings: We consider a 
ring with TV sites, see Fig. [21 For notational conve- 
nience we denote with d-f^rn n ~ 0, . . . , Af the fcrmionic 
operators in the upper arm. d^n, n = 0, . . . , N ~ M 



the fermionic operators in the lower arm, and identify 
d^o = d^o, d^M = dij,N-M- In terms of the operators 
d-fim rf^Tt the device Hamiltonian reads 

M-l M'-l 

Hdcv = t ^ (4„c?frn+i+h.c.)+t ^ (d[„dj;„+i+h.c.), 

ri=0 n=0 

(14) 

where M' = N — M and the hopping t is the same as 
in the leads. For simplicity we also set = tR = t 
in the tunneling Hamiltonian. We employ the TBWGA 
to calculate , e.g., the right-going eigenstates tAeL- Let 
i'shij) = \/ VL{£){e.'^^^ + i?eLe~''^^) be the amplitude on 
lead L {j < 0) and ipeL{j) = \fvL{e)T^Le^^i be the am- 
plitude on lead R (j > 0). Similarly, the wavefunction 
on the upper arm of the ring has the form ipgLln) = 
^i^i,(£)(A£ie*'=" + Beie-^'=") with Q<n< M while on 

the lower arm i'eL{n) = yJvL{e){CeLe'^" + L'^ie"'*'") 
with < n < N ~ M. According to the TBWGA the 
coefficients of ilJeL are solution of 



/I 





e 


e 



\kM 



ifc(M-l) 



ikM 
ik 

ik(M-l) 



Jk{N-M) 
Ak 



-ik(N-M) 
-ik 



Jk{N-M-l) ^-ik(N-M~l) Q 






-1 
-1 





ik 



ReL 



/I 
1 




\0 



(15) 



r 



Left-going states tp^R can be computed in a similar 
manner and it is straightforward to show that T^j^ = 
T^R = Tg, which is the transmittance of the system. The 
conductance through a bond is given by Eq. (|13p . For 
any bond in the leads the conductance is simply = 
iT'ePs'o where go = l/(27r) is the quantum of conductance 
for spinless electrons. The conductance on the upper and 
lower arm of the ring are G^^ = (I^^Lp — \Bi.L\'^)go and 
G^e = (|GeLp — |Z)eLp)go- From the above system of 
equations we obtain the following analytical solution for 
the transmittance and the local conductance 



jki2M-N) gjj^2(^) sin(fcf ) cos(fc(f - M)) 

n(k) 



(16) 



Gffe 

50 



sin''(fc) sin(fcf ) cos(fc(f - M)) siTi{k{N - M)) 



2|17(fc)|2 



Gjie _ sin^(fc)sin(fcf )cos(fc(f -M))sin(A:M) 



go 



(17) 



(18) 



with 



n{k) = 



,-2ik(N-M+2) 



16 

^2^kN _ 2g 



ikN / ^2ik 



(19) 



l)^ _ ^2ik{N~M+2) 



It is easy to verify that current conservation is fulfilled 
since G^£ -I- Gj^g = jTepgo = Gg. It is also worth empha- 
sizing that Ge is bounded between and go while G^^ 
and Gjie can be much larger than g^ and either positive 
or negative. We say that we are in a vortex regime if 
sign[G^e] = — sign[G4e]. From Eqs. (|17I18|) we conclude 
that for arms of different length a vortex regime always 
exists since 



G^e _ sin(fc(Ar - M)) 
G^ ~ sin(fcM) 



(20) 



Equation ((201) simple and transparent. The onset of a 
vortex occurs for those values of the incident momentum 
corresponding to an eigenenergy of either the isolated 
lower arm fc„ = nn/{N — M) or of the isolated upper 

arm k^m = wm/M. In the following we characterize the 
vortex regime and show how to predict the occurrence 
of large ring currents in coincidence with a vanishingly 
small current in the leads. 
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Vortex regime: As already pointed out in Refs. 
the transmittancc has two different fcinds of zeros. 
The numerator in Eq. (fTS)) vanishes either for energies 
that exactly match an eigenvalue of the isolated ring, i.e., 



(C) 

of G-ffe vanishes for k = while the numerator of G 



for kn 



(M) 



2m: /N (matching momenta) or for energies 



at which there is perfect destructive interference at the 

(X) 

right interface, i.e., for km = (2m + 1)tt/{N — 2M) {in- 
terference momenta). At these points the transmittancc 
vanishes provided the denominator Q,{k) ^ 0. From the 
linear system in Eq. (fT5|) one can easily show that ^{k) 
cannot be zero at the interference momenta since the 
wavefunction at the right interface vanishes and hence 
Tj vanishes as well. As the numerator of goes to 

(X'] 

zero as ~ [k — kn ) the denominator 0(fc) has to be 
finite. On the other hand, the denominator can van- 
ish at the matching momenta for special values of M/N. 

Expanding f2(fc) around fci'^'' one finds to first order 

n{k) « _( e4»A/»^/jV _ 1)2/16 _^ _ kiM)^^ ^j^j^ 

|7(n)| > ^iV(2iV - 3) + M(4iV + 2) > 0. Thus n{k) 

approaches zero as (fc - ki^^) for integer 2Mn/N. In 
these cases the simple zero of the denominator cancels 
the simple zero of the numerator and is finite (unless 
ki^^ is an interference momentum as well). Our condi- 
tion for the cancellation of zeros include the symmetric 
case M = N/2 already discussed in Ref. d as well as 
other cases, see Fig. [3] where = 16 and M = 6 and the 
zero at 7r/2 is cancelled by the denominator (the conduc- 
tance Ge has the same zeros as the transmittancc since 



00 
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FIG. 3: (Color online) Local conductance G^^/go (red), 
Giyjgo (blue) and Se = sign[sin(fc(A - M))/ sm{kM)] 
(black), see Eq. (|20p . versus the incident momentum k. The 
function Se is 1 in the laminar regime and -1 in the vortex 
regime. Vertical dashed lines are drawn in correspondence of 
k^k-^ (M), k = k^ (I), fc = fc^ (fr), fc = fc" (4), and of the 
zeros of fl{k) (D). 



From Eqs. (|17ll8p one can see that a zero of G^ implies 
a zero of both G-o-e and G^.^. Novel zeros, however, exist 
for the local conductances. Specifically, the numerator 



vanishes for k = km , 
G^e is zero at k 



see Fig. [31 Thus, the conductance 
k^f'^ unless k'if'^ is also a zero of the 
denominator. In this latter case G^^ remains finite since 

it exists an integer m such that fci^'' = fc^^ and hence 
sin(Arfc/2)sin(fc(Af - M)) - (fc - k^^^)^ which cancels 
the double zero of |f2(fc)pfi^ In Fig. [3] this cancellation 
takes place at = 7r/2. A similar reasoning apply to 
G^e. To summarize, the local conductances G^^, Gj^g as 
well as the transmittancc can have simple or double 
zeros while G^ double or quadruple zeros. 

In accordance with the above analysis the onset of the 
vortex phase occurs for incident momenta k = k^^^^ ^ 
k^'-^ or k = /c^-^) ^ k^^\ Below we show that in the 
vortex regime there are special values of the incident mo- 
mentum yielding a circular ring current much larger than 
the current in the leads. In order to quantify this effect 
we define the vortex function 



V{e) 



G 



G. 



4e 



G, 



(21) 



The modulus of V{£) is in the range (0,1) when G^e and 
G^e have the same sign (laminar regime) and is greater 
than one in the vortex regime. At the onset either G^e 
or G^e vanishes and |V^(e)| = 1- The vortex function 
diverges at the zeros of the total conductance and the 
magnitude of the vortex can be classified according to 
the nature of zeros of Gg. For double zeros of Gg (single 
zeros of T^) the difference G-i^e — G^e has a single zero 
since a double zero would imply k^^"^ ~ k^'-'^ which in turn 
implies rt{k) = 0. Expanding the conductances around 
the single zero of wc find V{e) ^ l/(£ — £s). The 
current fiowing in the ring changes direction as e crosses 
Eg ■ It is worth noting that in the neighborhood of Eg the 
derivative |G:J^^^^| can be very large (see Fig. i aX kp = 

Stt/S) and IG-j^/^^l ~ .go >> Gg ~ 10^ ^50- Even more 
striking is the behavior of V{£) around a double zero £d 
of Tj. In this case G^e — Gj^e also has a double zeroi^ 
and hence V{e) ^ l/(e — e^)^- In this case the vortex 
does not change sign as e crosses Ed- Both divergences 
arc rather remarkable since one could naively expect that 
[G^/^el < Ge in accordance with Kirchoff's current laws 
of classical electromagnetism. 

We also explored the vortex regime for nonzero onsite 
energies Edcv on the ring and different couplings iL/fl- 
Varying ti/jj alters the shape of G-i^/i^ but preserves both 
the position and the nature of the zeros. On the contrary, 
the position of the zeros changes as ejcv is varied but a 
vortex regime still exists. 

In conclusion wc have generalized the WGA to TB 
models and obtained a Landauer-like formula for the cur- 
rent density through a generic bond. The TBWGA re- 
quires the same computational effort as the WGA. Em- 
ploying the TBWGA in combination with the obtained 
expression of the local conductance we showed how to 
predict the occurrence of large vortex currents. The ex- 
istence of a vortex regime is rather robust and has to be 
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attributed to the nontrivial topology of the system rather Acknowledgment: S.B. acknowledges partial support 
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